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Constraining the Modified Newtonian Dynamics from spherically 
symmetrical hydrodynamic accretion 
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ABSTRACT 

The Modified Newtonian Dynamics (MOND) is an alternative to the dark matter assumption 
that can explain the observed flat rotation curve of galaxies. Here hydrodynamic accretion is 
considered to critically check the consistency and to constrain the physical interpretation of 
this theory. It is found that, in case of spherically symmetrical hydrodynamic accretion, the 
modified Euler's equation has real solution if the interpretation is assumed to be a modification 
of the law of dynamics. There is no real solution if it is assumed to be an acceleration scale 
dependent modification of the gravitational law. With the modified Euler's equation, the steady 
state mass accretion rate is found to change up to ~ 15%. The astrophysical and cosmological 
implications of these results are also discussed. 
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1 INTRODUCTION 



The rotation curve of spiral galaxies can not be explained in terms 
of gravitational potential of onl y the visible mass with any rea- 
sonable mass to light ratio (e.g. feubin & Fordlll970l Il98(t Isofuel 



19961 : ISofue & Rubinl feoOl : Spa no et alj|2008h . This indicates the 
presence of a significant fraction of gravitating matter with very 
high mass to light ratio. Various other observational and theoretical 
constraints from velocity dispersion of elliptical galaxies, baryon 
fraction in galaxy clusters, gravitational lensing, structure forma- 
tion, CMB_p_ow^r_^pec^rtmi, obs^n^ion_rf Lyman a forest etc . 
(e.g. iFaber & Jackson! [l976l IWu et al.lll998l: ISpringel et al.ll2005l; 



Clowe et al.1 120061 : iMassev et all 120071 : Iviel. Bolton & Haehneltl 
2009) also point out that a significant fraction of the mass of the 



Universe is in the form of the dark matter wi th very little or no elec- 
tromagnetic interaction (Be rtone et al.l2005l) . Though the dark mat- 
ter concept is widely accepted to be an explanation of these obser- 
vations, there is no general agreement on the composition and var- 



successful theory in this category is the theory of the Modified 
Newtonian Dynamics (MOND) propose d to explain rotat ion curves 
without and "hidden mass hypothesis' ' jMilgrom|[T983aHbT) . Essen- 
tially, the proposal of this theory is that the acceleration due to grav- 
itational force is not linearly proportional to the force at very low 
acceleration limit. This simple modification has so far successfully 
explained most o f the observed rotation curves for different types 
of galaxies (e.g. IScarpal 120031: ISwaters et al.1l2010h . and the few 
cases where it fails can be explained in terms of inadequate data, 
larg e asymmetries in the velocity fi eld and observational uncertain- 
ties dMilgromlfl 99 ll : IS watersll 1 9991) . The MOND pr edictions have 
been i ndependently tested against observations (e.g. iFamaev et al.l 
120071 : iTiret & Combesll2007l : ISanders & McGaughll2002l and ref- 
erences therein) and found to be consistent. Though there have 
been some criticism of this theory followed by the claim of a di- 
rect proof of dark matter at c luster scale from weak lensing ob- 
servations dciowe et al. | l2006h. there are ways o f accommodating 



Angus et al. l200d l2007h . and the issue is 



Vittorio&Silkll984 


; Davis et al.ll985l:IUmemura & Ikeuchil 19851: 


Navarro et al.lll996t 


de Blokll2005l). Though there have been indi- 



cations from some of the direct and indirect detection experiments 
of different types of dark matters, no conclusive results showing 
firm detection is still reported. 

However, over the time, a variety of alternative theories, with 
some modification of either laws of motion or law of gravita- 
tional force, have been suggested to explain some of the obser- 
vations without invoking the assumption of any dark matter (e.j 



Milgrom|[l983l:lBekenstein & Milgromlll 984 Isanderslll 9861: iFahn 1 
ISanderslll997l : iBrownstein & M offat 200cJ). A particularly 
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these observations (e.g. 
far from being settled. 

Whether or not MOND is an alternative to the dark matter 
scenario, there is no doubt that a complete theory of dark mat- 
ter must explain this success of MOND in predicting and explain- 
ing such observations. It is, hence, increasingly important to crit- 
ically verify its consistency at all scales. There are proposals to 
test predic tions of MOND or as s ociated theories in a variety of 
ways (e.g. Ide Lorenci e t al. 20091: iTrenkel et ail 1201 Oh . In this pa- 
per, I have considered the astrophysical case of spherical accretion 
to check the consistency and to constrain the physical interpretation 
of MOND. The background is outlined in Section [J2] and the out- 
come in MOND regime for different modifications are described 
in Section Sj3] Section Sj4]contains discussions on the results, and I 
present conclusions in Section Sj5] 
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2 BACKGROUND 

In MOND, the Newtonian equation of dynamics F = ma is mod- 
ified to F — mpa, where p = p(a) is a dimensionless parame- 
ter. This modification is significant at very low acceleration regime 
(below an acceleration of ao ~ 10 -10 m s _1 ). It is proposed that 
p ~ a I ao for a < ao and p — 1 for a > ao, and the exact form of 
p(a) may not have any serious consequences. Thus the motion due 
to gravitational force will be governed by 



GMm , 

- — r = mpa. 



(1) 



For the low acceleration at a large distance from the centre of a 
galaxy, this will im ply a = v 2 /r = yJGMao/r giving rise to a flat 
rotation curve. See Sande rs & McGaughl (2002) for a comprehen- 
sive critical review of the theory, a more generalized formulation 
and its i mplications, its ob servational supports and other details. 
Also see iBekensteinl d2004l) and references therein for the details 
of the "MOND inspired" relativistic, generalized theory of gravita- 
tion called TeVeS. For the purpose of this work, I will only use this 
theory to be phenomenological as summarized in equation ([T} to 
investigate its possible implications in the case of spherically sym- 
metrical hydrodynamic accretion. 

Note that in equation (TJ, p can be written as a modification 
of either the inertial term (F g = mpa) or modification of the grav- 
itation force term (F g /p = ma). Though they lead to the same re- 
sult for the rotation curve, the physical implication is significantly 
different. There may be situations where these two interpretations 
leads to drastically different results. Note that this phenomenologi- 
cal description of the dynamics is consist ent with the nonrel ativistic 
limit of TeVeS with spherical symmetry I IBekens tein 2004). 



3 ACCRETION IN MOND REGIME 

Here I have considered spherical accretion to check if, in the 
MOND regime, there is any change of physical conditions from 
that of the Newtonian scenario. In the Newtonian case, the gov- 
erning equations for spherically symmetrical hydrodynamic steady 
state accretion are the continuity equation and the Euler's equation 

i_J_ (pwr 2 )=0 (2 ) 



du 
dr 



ldP _ GM 
p dr r 2 



(3) 



where u(r) is the radial inward velocity, P(r) and p(r) are pres- 
sure and density related by an equation of state P = Kp 1 and 
polytropic index 7, and M is the mass of accretor. The sound speed 
in the medium c s (r) is related to P and p as c 2 = dPj dp. 

Starting with a boundary condition poo and c soo at a very 
large distance from the central accretor, equations (|2] and can 
be solved for a given M and 7 to derive the steady state den- 
sity and velocity profile. The solution of astrophysical interest is 
an unique transonic solution with a mass accretion rate of M = 
47rA(7)pooG 2 M 2 /Cs 00 , where A is a dimensionless constant. Note 
that the acceleration at the sonic point is of the order of ci^/GM. 
For a solar mass accretor and a typical ambient sound speed of 
~ 10 km s" 1 , this is ~ 10~ 7 , very much comparable to the MOND 
acceleration constant ao. Keeping this in mind, it is useful here to 
introduce a dimensionless parameter ao = ao/(ct,oa/GM) and to 
study the behaviour of the system for different values of this pa- 
rameter. 



In case of the hydrodynamic accretion, the dynamics is gov- 
erned by the interplay of three terms - pressure, gravitational force 
and inertia. In the aforementioned two different interpretations of 
the MOND modification, viz. modification of gravitational force 
and modification of the Newtonian dynamics, these three terms 
change in different ways. Thus, considering this case of hydrody- 
namic accretion gives us a chance to study any possible difference 
that may arise in these two interpretations. 



3.1 Modification of gravitational force 

In this case, we consider the modification of the gravitational force 
of the form Fg = F g /p = ma. As mentioned earlier, this can 
be derived from TeVeS at no nrelativistic limit assuming spherical 
symmetry dBeken stein 2004). With this modification, the continu- 
ity equation will remain unchanged from equation (0, but the Eu- 
ler's equation will be modified to 



du 
dr 



ldP 

p dr 



GM 
pr 2 



(4) 



where p = p(a/ao) and a = y/GMao/r. Note that, in MOND 
regime, with p(a/ao) = a/ao < 1, the gravitational term will be 
\/GMao/r, whereas in the low acceleration regime with mu — 
1, it will be same as the regular Newtonian term GM/r 2 . Using 
equation 10 and the equation of state, P and p can be eliminated 
from equation ((4) in terms of it, and can be re-written as 



(5) 



(u 2 -c 2 s )du _ 2ci GM 
u dr r pr 2 

With an initial condition pi and a^ at a small radius rj, equation 0, 
along with the equation of state and the continuity equation, can be 
solved numerically for the density and velocity profile. Note that 
for typical astrophysical condition, a at small radius is significantly 
larger than ao. So, for the inner region, the MONDian solution will 
not differ from the Newtonian solution. Here also, like the Newto- 
nian case, velocity at r ; should have a unique value for the solution 
to pass through the sonic point and to give the right accretion rate. 
For a given 7 and ao, equation can be solved with same initial 
conditions for both the Newtonian and the MONDian case. For the 
Newtonian case, at large radius, p and c s tends asymptotically to a 
constant value poo and c soc respectively, and the velocity u ~ r~ 2 
tends to zero. Interestingly, as shown in figure (0, for MONDian 
case, equation 10 does not have this asymptotic solution at large 
radius. Figure ([T} shows the Newtonian and MONDian solutions 
for 7 = 7/5 and 2o ~ 0.3, 1.0 and 3.0. The top and bottom panels 
show the density and the velocity field respectively. Both the den- 
sity and the velocity profile deviate from the Newtonian solution 
and diverges away from the asymptotic solution at large radius. In 
these plots, the density and the velocity values are scaled by poo 
and c aoo (of the Newtonian solution), and the radius is scaled by 
r c = GM/c 2 soo. 

One can understand this result analytically with the following 
arguments. In the MOND regime with p — a/ao < 1, equation (0 
can be written in terms of u = ^ as 



VGA-Jap 



(u 2 



(6) 



Since both the term in numerator scales as 1/r, with decreasing c 3 , 
u changes sign at large radius. Thus, u/c s is no more a monotoni- 
cally decreasing function with increasing radius, and the MONDian 
solution diverges from the Newtonian solution. It implies that there 
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Newtonian - 
MOND1 a n /(c% /GM) = 0.3 - 
MOND1 aJc' IGU| .1.0 
MOND1 a /(c 4 s /GM) = 3.0 



Newtonian 

MOND1 aJc' /GM) = 0.3 

MOND1 a n /(c 4 ,„/GM) = 1.0 

MOND1 a / c 4 s /GM = 3.0 



Newtonian 

MOND2 R1 - - 
MOND2 R2 



Newtonian 

MOND2 R1 — 
MOND2 R2 



Figure 1. Density and velocity profile in the MOND regime for the modified 
gravitational force case. Here 7 = 7/5 and So ~ 0.3, 1.0 and 3.0 (for c 300 
of the Newtonian solution). Solid line is for the Newtonian solution. 

is no solution of the flow with p, u and c s having the right physical 
boundary condition. 



Figure 2. Typical spherical accretion density and velocity profile in the 
MOND regime for the modified dynamics case. Here 7 = 7/5 and ao 1 
(for c S tx) of the Newtonian solution). Solid line is for the Newtonian solu- 
tion. Other two broken lines correspond to the two different roots of equa- 
tion (9j. One of the roots is imaginary at large radius but the other one 
satisfies physical boundary condition requirement. 



3.2 Modification of dynamics 

In this case also, the continuity equation will remain as in equation 
l|2}. But since F = mpa, the modified Euler's equation will be 



du 



ldP GM 

p dr r 2 



(7) 



Following a similar numerical analysis as in the earlier case, one 
can solve for u(r) and p(r) from equations ((2) and (|7J. As shown 
in figure ((2), one of the solutions in this case starts to deviate from 
that of the Newtonian one at a large radius, and have asymptotic 
value of poo and c soo lower than that of the Newtonian value with 
same mass accretion rate. The other solution becomes imaginary 
after a certain distance. Figure l(2) shows the Newtonian and these 
two MONDian solution for 7 = 7/5 and 3o«l. The deviation of 
Poo from that of the Newtonian one depends on the value of both 7 
and ao- 

Eliminating P and p using equation {2} and the equation of 
state, equation 10 can be rewritten as 



pu 



c, du 



'2c- 



dr 



GM 



(8) 



In the MOND regime, where p = a /ao = — ^ , this also be- 
comes a quadratic equation in u 



u 2 ,2 2cl GM 

— u H u + ( — ) = 

ao u r r z 

where the condition of a real solution for v! is 



a 2cj _ GM 

d ^ 4 r ~ r 2 ' 



(9) 



(10) 



This is an upper limit condition to u and does not contradict the 
required physical condition of velocity tending to zero at large ra- 
dius. Thus, in this case, there exist a physically meaningful solution 
where u tends to zero at large radius whereas density and sound 
speed asymptotically tend to p^ and c soo respectively. 

Effectively, for this interpretation, the general nature of the 
solution in MOND regime does not change from that of the Newto- 
nian solution. In details, however, the mass accretion rate changes 
to M — 4iv\(-f,ao)p<xiG 2 M 2 /c soo , where A is a dimensionless 
factor. This is further investigated by evaluating A for different 7 
and 5o. As shown in figure ((3}, the steady state mass accretion 
rate may change by up to ~ 15%. However, in case of astrophys- 
ical accretion, the accretor mass, ambient density and sound speed 
are often not so well determined to observationally distinguish this 
change between Newtonian and MONDian accretion. 
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Figure 3. Mass accretion rate parameter A as a function of acceleration 
parameter So for different equation of state with polytropic index 7. 



4 DISCUSSIONS 

The results derived in Section fj3]is for the spherically symmetri- 
cal and non-magnetized hydrodynamics accretion. While consider- 
ing these results, it should be kept in mind that, in reality, astro- 
physical accretion is a complex phenomenon. These assumptions 
of steady state, spherical symmetry and no importance of angu- 
lar momentum, and ignoring the possible role of self-gravity and 
magnetic field are for the simplicity of this semi-analytical investi- 
gation. However, in the low acceleration regime far away from the 
accretor, the effect of magnetic field, self gravity and angular mo- 
mentum may be negligible, and the large scale accretion may be 
approximated as hydrodynamic, non-magnetized, spherical accre- 
tion onto a central compact accretor. Thus, the exact solution may 
get modified due to these complications, but the general nature of 
the solution is expected to remain unchanged. 

The other point to note is that the pressure term in the Euler's 
equation is assumed to be not affected by the MONDian modifica- 
tion. Fundamentally, this term, which arises form the random mo- 
tion and change of momentum of the particles, is expected to get 
modified in a similar way as the inertia term. In that case, equation 
will be reduced to equation 10, and, there will be no real solu- 
tion in the MOND regime. The way out is to consider equation ([T} 
as the governing equation not for random motion but only for bulk 
motion of system with symmetry. This scenario is consistent with 
the observation that the random motion and acceleration of the gas 
does not alter the galaxy rotation curve either. 



5 CONCLUSIONS 

For the case of spherically symmetrical hydrodynamic accretion in 
MOND regime, it is shown here that physically meaningful solu- 
tion exists only for the interpretation of the modification of dynam- 
ics but not of the gravitational law. At a phenomenological level, 
this modification should be not for random motion but for bulk mo- 
tion only. Given the uncertainty on various parameters, the change 
of accretion rate is not significant to distinguish between the New- 
tonian and the MONDian scenario. 
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